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Abstract 

The semi- inclusive difference asymmetry A'^'^~^ for hadrons of opposite charge has been 
measured by the COMPASS experiment at CERN. The data were collected in the years 
2002-2004 using a 160 GeV polarised muon beam scattered off a large polarised ^LiD target 
in the kinematic range 0.006 < x < 0.7 and 1 < < 100 (GeV/c)^. In leading order 
QCD (LO) the deuteron asymmetry measures the valence quark polarisation and 

provides an evaluation of the first moment of Au„ + Ady which is found to be equal to 
0.40 ± 0.07 (stat.) ± 0.06 (syst.) over the measured range of a; at ^ 10 (GeV/c)^ When 
combined with the first moment of gf previously measured on the same data, this result 
favours a non-symmetric polarisation of light quarks Au — — Ad at a confidence level of two 
standard deviations, in contrast to the often assumed symmetric scenario Au = Ad = As = 
As. 
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The COMPASS experiment at CERN has recently pubhshed an evaluation of the deuteron 
spin-dependent structure function gfi^x) in the DIS region, based on measurements of the spin 
asymmetries observed in the scattering of 160 GeV longitudinally polarised muons on a longi- 
tudinally polarised ^LiD target [Ij. These measurements provide an accurate evaluation of the 
first moment of gi for the average nucleon in an isoscalar target = {gf + g^)/2 

rf (Q2 = io (GeV/c)2) = / g^{x,Q^ = 10 {Ge\/cf)dx = 0.051 ± 0.003 (stat.) ± 0.006 (syst.) 

^0 

(1) 

from which the first moment of the strange quark distribution can be extracted if the value of 
the octet matrix element (as = 3F — D) is taken from semi-leptonic hyperon decays^]] At LO 
in QCD the strange quark polarisation is given by 

As + As = Srf - = -0.09 ± 0.01 (stat.) ± 0.02 (syst.) (2) 
at Q2 = 10 (GeV/c)2. 

Since quarks and antiquarks of the same fiavour equally contribute to 171, inclusive data 
do not allow to separate valence and sea contributions to the nucleon spin. We present here 
additional information on the contribution of the nucleon constituents to its spin based on 
semi-inclusive spin asymmetries measured on the same data as those used in Ref. [1]. 

The semi-inclusive spin asymmetries for positive and negative hadrons h'^ and h~ are 
defined by 

h+ _ h+ h- _ h- 

Ah+ _ '^Ti ^TT Ah- _ ^Ti ^TT /o^ 

'^n "^"TT Ti TT 

where the arrows indicate the relative beam and target spin orientations. 

The data used in the present analysis were collected by the COMPASS collaboration 
at CERN during the years 2002-2004. The event selection requires a reconstructed interaction 
vertex defined by the incoming and scattered muons and located inside one of the two target cells 
[2]. The energy of the beam muon is required to be in the interval 140 GeV < < 180 GeV and 
its extrapolated trajectory is required to cross entirely the two cells in order to equalise the fiuxes 
seen by each of them. DIS events are selected by cuts on the photon virtuality (Q^ > 1 (GeV/c)^) 
and on the fractional energy of the virtual photon (0.1 < y < 0.9). Final state muons are 
identified by signals collected behind the hadron absorbers. The hadrons used in the analysis are 
required to originate from the interaction vertex and to be produced in the current fragmentation 
region. The latter requirement is satisfied by selecting hadrons with fractional energy z > 0.2. 
In addition an upper limit z < 0.85 is imposed in order to suppress hadrons from exclusive 
diffractive processes and to avoid contamination from muons close to the beam axis which escape 
identification by the muon filters. The hadron identification provided by the RICH detector is 
not used in the present analysis. The resulting sample contains 30 and 25 million of positive and 
negative hadrons, respectively. 

The target spins are reversed at regular intervals of 8 hours during the data taking. The 
spin asymmetries are obtained from the numbers of hadrons collected from each target cell during 
consecutive periods before and after reversal of the target spins, following the same procedure as 
for inclusive asymmetries [3]. They are listed in Table [Hand also shown in Fig. [T] as a function of 
X, in comparison to the SMC [11[5] and HERMES [6j results. The consistency of the results from 
the three experiments illustrates the weak dependence of the semi-inclusive asymmetries. 
The COMPASS results show a large gain in statistical precision with respect to SMC, especially 
in the low x region (x < 0.04), while at larger x the COMPASS errors are comparable to those 
of HERMES. The systematic errors, shown by the bands at the bottom of the figure, result from 



At the precision of the experiment the value of is unchanged when the evolution of the measured values 
gi{xi, Ol) to a common is done at LO or at NLO in QCD. 




Figure 1: Hadron asymmetries Ay' (left) and (right) measured by COMPASS, SMC [5] and 
HERMES [6j experiments. The bands at the bottom of the figures show the systematic errors 
of the COMPASS measurements. 



different sources. The uncertainty on the various factors entering in the asymmetry calculation 
(beam and target polarisation, depolarisation factor and dilution factor) leads to a relative 
error of 8% on the asymmetry when combined in quadrature. The uncertainty due to radiative 
corrections is smaller than in the inclusive case due to the selection of hadronic events and does 
not exceed 10~^ in any x bin. The presence of possible false asymmetries due to time-dependent 
apparatus effects has been studied in the same way as for the inclusive asymmetries: the data 
sample has been divided into a large number of subsamples, each of them collected in a small 
time interval. The observed dispersion of the asymmetries obtained for these subsamples has 
been found compatible with the value expected from their statistical error. This allows to set an 
upper limit for this type of false asymmetries at about half of the statistical error. Asymmetries, 
obtained with different settings of the microwave frequency used for dynamic nuclear polarisation 
of the target, have also been compared and did not reveal any systematic difference. 

Under the common assumption that hadrons in the current fragmentation region are pro- 
duced by independent quark fragmentation, the semi-inclusive asymmetries A^^ can be 
written in LO approximation as 



^ , _ E,egAg(x,Q^)Z),^(.,Q^) 



where Aq{x,Q ) and q{x,Q ) are the polarised and unpolarised parton distribution functions 
(pdf's) and D^{z,Q'^) the fragmentation function of a parton q into a hadron h. The above 
formula does not account for the full complexity of the hadronisation process as described, 
for instance in the Lund string fragmentation model [7], and its validity in low energy fixed 
target experiments has been questioned [8]. It has nevertheless been shown to hold as a good 
approximation at the energy of COMPASS [9]. 

In addition to purely experimental effects such as (x, z) correlations in the spectrometer 
acceptance, a z dependence of the semi-inclusive asymmetries A^~^,A^~ could reveal a break- 
down of the independent fragmentation formula (Eq. In order to check the possible presence 
of this effect, we have re-evaluated the asymmetries for each x bin subdivided into three inter- 
vals of z. Within their statistical precision the obtained values do not indicate any systematic z 
dependence. 

In the present analysis we use the "difference asymmetry" which is defined as the spin 
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Table 1: Values of , A^ and A^ with their statistical and systematical errors as a 
function of x with the corresponding average value of Q^. 



asymmetry for the difference of the cross sections for positive and negative hadrons 

.h+-h- ^ ^^n ' '-'^TT '^TT ) 

The difference asymmetry approach for the extraction of helicity distributions, introduced in 
Ref. [lO], has been used in the SMC analysis ^ and been further discussed in [11^ I12j . In LO 
QCD and under the assumption of isospin and charge conjugation symmetries, the fragmentation 
functions cancel out from A" . In addition, in the case of an isoscalar target and assuming 
As = As, the difference asymmetries for pions and kaons are both equal to the valence quark 
polarisation 

Al'--=Ar--- = ^^^^i±4^, (6) 



where we introduce the valence quark distributions Qy = q — q- Since kaons contribute to the 
asymmetry in the same way as pions, their identification is not needed, allowing to reduce the 
statistical errors. The difference asymmetry for (anti)protons A^^^ has also the same value but 
under more restrictive assumptions and is more likely to be affected by target remnants. Since 
protons and antiprotons account only for about 10% of the selected hadron sample, the relation 



^h+-h- _ An^ + Ady 

Uy ± dy 



is expected to hold as a good approximation in the present analysis. Monte Carlo studies using 
POLDIS [13j and Lund string hadronisation show that the asymmetries A'^^{x) closely follow 
the trend of A'^~'^ (x) with a difference never exceeding 0.02. In addition the semi-inclusive 
asymmetries A^~^ are found to be very close to the expected values (Au^ ± A(i^)/(u^ ± d^) 
defined by the input parametrisations in the Monte Carlo simulation with the largest difference 
(< 0.05) appearing in the two highest intervals of x. 

At higher order in QCD the difference asymmetries still determine the valence quark 
polarisation without any assumption on the sea and gluon densities jllj . Fragmentation functions 
no longer cancel out but their effect is expected to be small [12]. 

The relation between the difference asymmetries of Eq. ([5]) and the single hadron asym- 
metries of Eq. ([3]) is 

■4'*-'- = r^C^'* - rA^-), with . = = (8) 

1 r +a^^ t 



The ratio of cross sections for negative and positive hadrons r depends on the event kinematics 
and is obtained as the product of the corresponding ratio of the number of observed hadrons 
N~ /N~^ by the ratio of the geometrical acceptances a'^ /a~ 
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h- 



N- 



a ' 
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(9) 



The ratio of the number of negative to positive hadrons (Fig. [2l left) decreases with increasing x. 
This ratio is subject to acceptance corrections because positive and negative hadrons, produced 
at the same angle, cross different regions of the spectrometer. To this end LEPTO generated 
Monte Carlo events have been processed through the program simulating the COMPASS spec- 
trometer performance [2j and reconstructed in the same way as the data. The acceptances a"*" 
and a~ are indeed found to be different. The ratio a~ /a^ , which is about 1.0 at low x, increases 
for X > 0.1 reaching ~1.12 in the highest x bin. Bin migration was found to be negligible. The 
corrected cross section ratio /a^^ is also shown in Fig. [2j 
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Figure 2: Left: The ratio /a^* before (triangles) and after acceptance corrections (circles). 
Right: The difference asymmetry, A^~^ , for unidentified hadrons of opposite charges, as a 
function of x at the of each measured point. 



The resulting values of the difference asymmetry as a function of x are shown in 

Fig. [2] (right) and listed with their statistical and systematic errors in Table [TJ The statistical 

which is approximately 0.20 over the measured range of x, 



h+-h- 



As can be seen from Eq. 



correlation between and 

is taken into account in the evaluation of the error of A 
a singularity appears when the cross section ratio becomes close to one, leading to infinite 
statistical errors. For this reason, we discard the lowest x bin used in the inclusive gi analysis [T] 
and take x = 0.006 as lower limit for the present analysis. The increase of A^^ for X > 0.1 
illustrates the increasing polarisation of valence quarks carrying a larger fraction of the nucleon 
momentum. 

Target remnants may affect current quark fragmentation at low values of the total hadronic 
energy W. In order to check the possible presence of such effects in our data, we have re-evaluated 
the difference asymmetries A'^^~^ (x) with the cut W > 7 GeV. The comparison of the obtained 
values with those quoted in Table 1 shows that the W cut only affects the two highest intervals 
of X where {x) is reduced by about 0.3 cistat- It will be shown below that these two 

values can be replaced by more accurate estimations so that the observed W dependence does 
not affect any further result. 

The polarised valence quark distribution Au^ + Ad^ is obtained by multiplying A'^ ~^ 
by the unpolarised valence distribution of MRST04 at LO - Here two corrections are applied, 
one accounting for the fact that although R{x,Q^) = at LO, the unpolarised pdf's originate 



from F2^s in which R = ctl/ctt was different from zero [15], the other one accounting for deuteron 
D-state contribution {ujd = 0.05 ± 0.01 |16]): 



{Uy + rft,)MRST 

^l + i2)(l - 1.5l^d) 



A' 



h+-h- 



(10) 



The LO parameterisation of the DNS fit [T7] has been used to evolve all values of Au^ + Ady to 
a common fixed at Qq = 10 (GeV/c)^ assuming that the difference between An„ + Ad^ at 
the current and at Qq is the same for the data as for the fit |Tj. The DNS analysis includes 
all DIS gi data prior to COMPASS, the partial COMPASS data on gi from Ref. [3] as well as 
the SIDIS data from SMC |^ and HERMES [6j. Two parameterisations of polarised pdf's are 
provided at LO, corresponding to two different choices of fragmentation functions, KRE [12] and 
KKP [20]. We have checked that the x dependence of the ratio o"'* /a^'^ (Fig. [2|) is fairly well 



reproduced by the LO MRST04 pdf's and the KKP fragmentation functions whereas the KRE 
parameterisation leads to a much weaker x dependence. For this reason we choose the fit with 
the KKP parameterisation. The largest corrections to x[Au„(x) + Ady{x)] are at large x and 
and do not exceed 0.03. The use of different fits (NLO fit of Ref. p7] or [18]) leads practically 
to the same results. The resulting values are shown in Fig. [3] (left). The DNS fit, also shown in 
the figure, is basically defined by the SMC and HERMES semi-inclusive asymmetries. Its good 
agreement with the COMPASS values {x^ = 7.7 for 11 data points) illustrates the consistency 
between the three experiments. 

The sea contribution to the unpolarised structure function F2 decreases rapidly with in- 
creasing X and becomes smaller than 0.1 for x > 0.3. Due to the positivity conditions |Ag| < q 
and \Aq\ < q, the polarised sea contribution to the nucleon spin also becomes negligible in this 
region. In view of this, the evaluation of the valence spin distribution of Eq. (jlOp can be replaced 
by a more accurate one obtained from inclusive interactions. Indeed at LO one obtains 



Auy + Adv 



36 



9i 



5 (1 - 1.5wd) 



2(An + Ad) + -{As + As) 
5 



(11) 



The values obtained by taking only the first term on the r.h.s. for x > 0.3 are also shown 
in Fig. [3l They agree very well with the DNS curve, which is based on previous experiments 
where the same procedure had been applied [3 [6]. The upper limit of the neglected sea quark 
contribution, derived from the saturation of the positivity constraint \Aq\ < q is included in the 
systematic error. 

The first moment of the polarised valence distribution, truncated to the measured range 

of x, 

fO.7 



Ft) {xj- 



[Auy{x) + Adtj{x)] dx, 



(12) 



derived from the difference asymmetry for x < 0.3 and from gf for 0.3 < x < 0.7, is shown 
in Fig. [3] (right). Practically no dependence on the lower limit is observed for Xmin < 0.03. We 
obtain for the full measured range of x 



F4O.OO6 < X < 0.7) = 0.40 ± 0.07 (stat.) ± 0.06 (syst.) (13) 

at Q2 = 10 (GeV/c)^ with contributions of 0.26 ± 0.07 and 0.14 ± 0.01 for x < 0.3 and x > 0.3, 
respectively. The uncertainty due to the unpolarised valence quark distributions (~ 0.04) has 
been estimated by comparing different LO parametrisations and been included in the systematic 
error. It should be noted that removing the factor (1 -|- R) in Eq. (10) would increase the 
value of F„ to 0.42 ± 0.08 ± 0.06. Our value of F„ confirms the HERMES result obtained at 
= 2.5 (GeV/c)^ over a smaller range of x and is also consistent with the SMC result which 
has three times larger errors (Table [2|). The factor (1 + R) was also used in the analyses of the 
previous experiments. 





x-range 


(GeV/c)2 


An^ + Ad^ 


Au± Ad 


Exp. Value 


DNS 


Exp. Value 


DNS 


SMC 


0.003 - 0.7 


10 


0.26 ±0.21 ±0.11 


0.386 


0.02 ±0.08 ±0.06 


-0.009 


HERMES 


0.023 - 0.6 


2.5 


0.43 ± 0.07 ± 0.06 


0.363 


-0.06 ±0.04 ±0.03 


-0.005 


COMPASS 


0.006 - 0.7 
- 1 


10 


0.40 ± 0.07 ± 0.06 
0.41 ± 0.07 ± 0.06 


0.385 


0.0 ±0.04 ±0.03 


-0.007 



Table 2: Estimates of the first moments An^±A(ii, and Aii±Ad from the SMC [5j, HERMES [6j, 
COMPASS data and also from the DNS fit at LO [17J truncated to the range of each experiment 
(lines 1-3). The SMC results were obtained with the assumption of a SU{3)f symmetric sea: 
An = Ad = As. The last line shows the COMPASS results for the full range of x. 



The difference between our measured value of r„ (0.006 < x < 0.3) and the integral of gi 
over the same range of x gives a global measurement of the polarised sea. Indeed, re-ordering 
Eq. (11) we obtain 



0.30 _ 

' dx = -0.02 ± 0.03 (stat.) ± 0.02 (syst.), (14) 



0.006 



(Au ± Ad) ± - (As ± As) 
5 



where the correlation between inclusive and semi-inclusive asymmetries has been taken into 
account in the statistical error. This result is compatible with zero but also consistent with the 
strange quark contribution of Eq. ([2]) and a vanishing contribution from the first term. It should 
be kept in mind that moments of sea quarks evaluated at LO have to be taken with caution 
because their values are small and thus comparable to the NLO corrections. 

The unmeasured contribution to for x > 0.7 estimated from the LO DNS parameterisa- 
tion of Ref. [T7] is 0.004 at = 10 (GeV/c)^. Its upper limit corresponding to the assumption 
A^^~^ = 1 for x > 0.7 is 0.007 according to the MRST04 parameterisation. 

The unmeasured low x contribution to Ty is expected to be negligible since the integral 
shows no significant variation when its lower limit is varied between 0.006 and 0.02. We thus 
estimate the first moment as 



r^(0 < X < 1) = 0.41 ± 0.07 (stat.) ± 0.06 (syst.). (15) 

The assumption of a fully flavour symmetric sea Am = Ad = As = As obviously leads to 
r^(0 < X < 1) = as. As shown in Fig. [3] (right), our experimental value is two standard 
deviations below the value of ag = 0.58 ± 0.03 derived from hyperon /? decays |21j . It has been 
suggested that a value of the valence contribution F^ smaller than as (as expected from the 
constituent quark models) could be a hint that a so far unmeasured part of the nucleon's spin 
resides at x = [22] . 

An estimate of the light sea quark contribution to the nucleon spin can be obtained by 
combining the values of F„ (Eq. (13)), F^ (Eq. (1)) and as 

Au + Ad = 3Ff - ^F, ± ^as (16) 

and the result is found to be zero (Table [2]) . Possible deviations from the nominal value of ag 
due to SU{3) f symmetry violation in hyperon decays are generally assumed to be of the order of 
10% [23| and are included in the systematic error. The zero value of An± Ad is in contrast with 
the non-zero value obtained for As ± As (Eq. ([2])) and suggests that Au and Ad, if different 
from zero, must be of opposite sign. Previous estimates by SMC and HERMES, also given in 
Table [51 are compatible with this hypothesis. The DNS parameterisation flnds a positive Au 
and a negative Ad, about equal in absolute value. Opposite signs of Au and Ad are predicted in 
several models, e.g. in Ref. [24j (see also [25] and references therein). Forthcoming COMPASS 
data on a proton target will provide separate determinations of Au and Ad. 




Figure 3: Left: Polarised valence quark distribution x[Auv{x) + Adv{x)) evolved to = 
10 (GeV/c)^ according to the DNS fit at LO [ITj (line). Three additional points at high x 
are obtained from gf [T]. The two shaded bands show the systematic errors for the two sets of 
values. Right: The integral of Auy{x)+Ad^{x) over the range 0.006 < x < 0.7 as the function 
of the low X limit, evaluated at = 10 (GeV/c)^. 

In conclusion, we have determined at LO QCD the polarised valence quark distribution 
from the difference asymmetry for oppositely charged hadrons in DIS of muons on a polarised 
isoscalar target. Its first moment at = 10 (GeV/c)^ over the measured range of x (0.006-0.7) 
is found to be 0.40 ±0.07 (stat.) ±0.06 (syst.). This value disfavours the assumption of a flavour 
symmetric polarised sea at a confidence level of two standard deviations and suggests that Au 
and Ad are of opposite sign. 
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